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EVERY INTEGER GREATER THAN 454 IS 
THE SUM OF AT MOST SEVEN POSITIVE CUBES 

SAMIR SIKSEK 


Abstract. A long-standing conjecture states that every positive integer other than 
15, 22, 23, 50, 114, 167, 175, 186, 212, 

231, 238, 239, 303, 364, 420, 428, 454 

is a sum of at most seven positive cubes. This was first observed by Jacobi in 1851 on the 
basis of extensive calculations performed by the famous computationalist Zacharias Dase. We 
complete the proof of this conjecture, building on previous work of Linnik, Watson, McCurley, 
Ramare, Boklan, Elkies, and many others. 


1. Historical Introduction 

In 1770, Edward Waring stated in his Meditationes Algebraicae, 

Omnis integer numerus vel est cubus, vel e duobus, tribus, 4 , 5, 6, 7, 8, vel novem 
cubis compositus, ... 

Waring’s assertion, can be concisely reformulated as: every positive integer is the sum of nine 
non-negative cubes. Henceforth, by a cube we shall mean a non-negative cube. In the 19th 
century, numerical experimentation led to refinements of Waring’s assertion for sums of cubes. As 
noted by Dickson (1927), “At the request of Jacobi, the famous computer Dase constructed a table 
showing the least number of positive cubes whose sum is any p < 12000 In an influential Crclle 
paper, Jacobi (1851) made a series of observations based on Dase’s table: every positive integer 
other than 23 and 239 is the sum of eight cubes, every integer > 454 is the sum of seven cubes, and 
every integer > 8042 is the sum of six cubes. Jacobi believed that every sufficiently large integer is 
the sum of five cubes, whilst recognizing that the cut-off point must be far beyond Dase’s table, 
and he wondered if the same is true for sums of four cubes. He noted that integers e4, 5 (mod 9) 
cannot be sums of three cubes. Later computations by Romani (1982) convincingly suggest that 
every integer > 1290 740 is the sum of five cubes, and by Deshouillers, Hennecart, and Landreau 
(2000) that every integer > 7373170 279 850 is the sum of four cubes. 

Progress towards proving these observations of Waring, Jacobi and others has been exceedingly 
slow. Maillet (1895) showed that twenty-one cubes are enough to represent every positive integer. 
At the heart of Maillet’s proof is an idea crucial to virtually all future developments; the identity 
(r + a;) 3 + (r - x) 3 = 2r 3 + 6rx 2 allows one to reformulate the problem of representing an integer 
as the sum of a (certain number of) cubes in terms of representing a related integer as the sum of 
(a smaller number of) squares. Exploiting this idea, Wieferich (1908) proved Waring’s assertion 
(Wieferich’s proof had a mistake that was corrected by Kempner (1912)). In fact, the theoretical 
part of Wieferich’s proof showed that all integers exceeding 2.25 x 10 9 are sums of nine cubes. 
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Completing the proof required appealing to a table of von Sterneck (1903) (who extended Dase’s 
table to 40000), and applying what is now known as the greedy algorithm to reach the bound. 

Soon thereafter, Landau (1911) showed that every sufficiently large integer is the sum of eight 
cubes. This was made effective by Baer (1913), who showed that every integer > 14.1 x 233 6 « 
2.26 x 10 15 is the sum of eight cubes. Dickson (1939) completed the proof of Jacobi’s observation 
that all positive integers other than 23 and 239 are sums of eight cubes. Remarkably, Dickson’s 
proof relied on extending von Sterneck’s table to 123 000 (with the help of his assistant, Miss 
Evelyn Garbe) and then applying the greedy algorithm to reach Baer’s bound. 

In 1943 Linnik showed that every sufficiently large integer is the sum of seven cubes. A sub¬ 
stantially simpler proof (though still ineffective) was given by Watson (1951). Linnik’s seven 
cubes theorem was first made effective by McCurley (1984), who showed that it is true for inte¬ 
gers > exp(exp(13.94)). Ramare improved this to exp(205000) in 2005 and finally to exp(524) & 
3.72 x 10 227 in 2007. This bound is way beyond computer searches combined with the greedy 
algorithm. In Desliouillcrs et al. (2000), it is shown that every integer between 1290 741 and 10 16 
is a sum of five cubes. As observed by Ramare (2007), combining this with the greedy algorithm 
(Bertault et al., 1999, Lemma 3), we can easily deduce that every integer 455 < N < exp(78.7) « 
1.51 x 10 34 is the sum of seven cubes. 

There has been a number of partial results concerning sums of seven cubes. Bertault et al. 
(1999) show that every non-negative integer which is a cubic residue modulo 9 and an invertible 
cubic residue modulo 37 is a sum of 7 cubes. Boklan and Elkies (2009) show that every multiple 
of 4 greater than 454 is the sum of seven cubes, whilst Elkies (2010) shows the same for integers 
= 2 (mod 4). 

In this paper we complete the proof of Jacobi’s seven cubes conjecture, building on the afore¬ 
mentioned great works. 

Theorem 1. Every positive integer other than 

15, 22, 23, 50, 114, 167, 175, 186, 212, 231, 238, 239, 303, 364, 420, 428, 454 
is the sum of seven cubes. 

The programs that accompany this paper are available from http://tinyurl.com/zlaeweo. 
It is a pleasure to thank Alex Bartel, Tim Browning, John Cremona, Roger Heath-Brown and 
Trevor Wooley for stimulating discussions. 

2. The Main Criterion 

Let JE = exp(524) and JT' = exp(78.7). By the results of Ramare (2007) and of Deshouillers et al. 
(2000), it is sufficient to prove that every integer W < N < Jff is the sum of seven cubes. The 
results of Boklan and Elkies (2009) and Elkies (2010) allow us to restrict ourselves to odd integers 
N (our method can certainly be adapted to deal with even integers, but restricting ourselves to 
odd integers brings coherence to our exposition). In this section we give a criterion (Proposi¬ 
tion 2.2) for all odd integers N in a range K\ < N < K 2 to be sums of seven cubes. Most of 
the remainder of the paper is devoted to showing that this criterion holds for each of the ranges 
(9/10) n+1 J(h < N < (9/10)"J^ with 0 < n < 4226. This will complete the proof of Theorem 1 as 
(9/10) 4227 JT « 1.42 x 10 34 and Jf' » 1.51 x 10 34 . 

Theorem 2 (Gauss, Legendre). Let k > 0 be an even integer. There exist integers x, y, z such 
that 

( 1 ) 


x 2 + x + y 2 + y + z 2 + z = k. 
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Proof. Dividing by 2 we see that this is in fact the famous theorem, due to Gauss, that every 
non-negative integer is the sum of three triangular numbers. Alternatively, we can rewrite (1) as 

(2) (2x + l ) 2 + (2 y + l ) 2 + (2z + l ) 2 = 4 k + 3. 

As k is even, 4fc + 3 = 3 (mod 8 ); by a theorem of Legendre, every positive integer = 3 (mod 8 ) is 
the sum of three odd squares. □ 


Throughout this section m will denote a positive integer satisfying the conditions 

(i) m is a squarefree, 

(ii) 3 | m, 

(iii) every prime divisor of m/3 is = 5 (mod 6). 

Observe that m = 3 (mod 6 ). Moreover, for any integer TV, there is a unique integer t e [0, m) such 
that TV = 8 t 3 (mod m). Our starting point is a modified version of Lemma 3 of Watson (1951). 

Lemma 2.1. Let 0 < K i < K 2 be real numbers. Let m be a positive integer satisfying (i)-(iii) 
above. Let £ m , 5 m be real numbers satisfying 

(iv) 0 < £ m < S m < 1 , 

(v) K\ > (8 Sfn + 1/36 )m 3 + 3m/4, 

(vi) K 2 < ( 8 c^ + l/18)?n 3 + m/ 2 . 

Let K 1 < TV < K 2 be an odd integer. Suppose TV = 8 t 3 (mod m ) with t e [e m ■ ?n,S m ■ m ). Then TV 
is the sum of seven non-negative cubes. 


Proof. Write m = 6 r + 3. Let 

(3) k= N ~ 8t3 ~(r 2 +r+l). 

m 

The ciuantity k is an integer as TV = 8 1 3 (mod m), and even as (TV- 8 f 3 )/m and r 2 + r+ 1 are both 
odd. Observe that 


TV - 8 5^ • to 3 , 2 . , „ 

k > - (r +r+ 1) as t < S m - m 


m 

K\ - 8 Si ■ m 3 


(r 2 + r + 1) as TV > K\ 


K\ - 86 3 • m 3 m 2 3 


>0 


m 

by (v). 


- 7 ^- - — substituting r = (m — 3 )/6 


As k is non-negative and even, by the Gauss-Legendre theorem, there exist integers x, y, z 
satisfying (1). We shall make use of the identity 


(4) (r + 1 + x) a + (?’ - x) 3 + (r + 1 + y ) 3 + (r - y) 3 + (r + 1 + 2) 3 + (r - z) 3 - 

(6 r + 3 )(?’ 2 + r + l + x 2 +x + y 2 + y + z 2 + z). 

From the definition of k in (3) and the fact that m = 6 r + 3, we see that TV - 8 f 3 is equal to the 
right-hand side of the identity (4). Hence 

TV = (r + 1 + x) 3 + (r - x) 3 + (r + 1 + y ) 3 + (r - y ) 3 + (r + 1 + z) 3 + (r - z) 3 + (2 1) 3 . 

To complete the proof it is enough to show that these cubes are non-negative, or equivalently that 

-r - 1 < x, j/, 2 < r. 

This is equivalent to showing that 

-( 2 r + 1 ) < 2x + 1 , 2y + 1 , 2z + 1 < 2 r + 1 . 
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Now (2 y+ l) 2 , (2 z+ l) 2 > 1 and so from (2), we have (2x + l) 2 < Ak+ 1. It is therefore enough to 
show that 4/c + 1 < (2r + l) 2 or equivalently that k < r 2 + r. The following inequalities complete 
the proof: 

N — 8 

k - r 2 - r =-(2?’ 2 + 2 r + 1) from (3) 


N - 8e 3 • m 3 


- (2 r 2 + 2r + 1) as t > e m ■ m 


(2r 2 + 2r + 1) as N < K 2 


K 2 - • m 3 m 2 1 

m 18 2 

<0 by (vi). 


substituting r = (m — 3)/6 


This simple-minded lemma has one serious flaw. The condition K\ < K 2 together with condi¬ 
tions (iv), (v), (vi), imply 

<t<<4 + l/288-l/(32m 2 ). 

In particular, this forces the interval [e-m, 8-m) to have length < m / \/288 » 0.15m. On the other 
hand, the integer t appearing in the lemma (which is the cube-root of IV /8 modulo in) can be any 
integer in the interval [0,?n). Thus the lennna only treats a small fraction of the odd integers 
K\ < N < K 2 . Our key innovation over the works mentioned in the introduction is to use not just 
one value of m, but many of them simultaneously. Each value of to will give some information 
about those odd integers Ki< N < K 2 that cannot be expressed as sums of seven cubes; collecting 
this information will allow us to deduce a contradiction. 

Let 1 be a real number and m be a positive integer. Define the quotient and remainder 
obtained on dividing x by m as 

Q(x, m) = lyc/jnJ, R(:r, m) = x - Q(x, m) ■ m. 

In particular, R(;c,??t.) belongs to the half-open interval [0,m). If x e Z then R (x,m) is the usual 
remainder on dividing by m, and x = R (x,m) (mod m). Let e and S be real numbers satisfying 
0 < e < 5 < 1. Define 

(5) Bad(m,£,<5) = ■{ x e M : R(a;, to) e [0,m) \ \e ■ m, 6 ■ to) > = (J fcm+([0,m) \ [e-m, 5 -to)) . 

k =—00 

The reader will observe, in Lemma 2.1, if N is not the sum of seven cubes, then t e Bad(m,£, 5), 
which explains our choice of the epiphet ‘bad’. Given a set of positive integers W, and sequences 
£ = (£m)meWi § = (<f m )m£W of real numbers satisfying 0 < £ m < S m < 1 for all m e W, we define 

( 6 ) Bad(W,£,5)= H Bad (m,e m ,S m ). 

meVV 

To make the notation less cumbersome, we usually regard the values e m and S m as implicit, and 
write Bad(m) for Bad(? 7 i,£,d), and Bad(W’) for Bad(>V, £, 5). 

Proposition 2.2. Let 0 < K\ < K 2 be real numbers. Let W be a non-empty finite set of integers 
such that every element m e W satisfies conditions (i)-(iii). Suppose moreover, that for each 
to e W, there are real numbers e m , 5 m satisfying conditions (iv)-(vi). Let M = lcm(W). Let 
6 c [ 0 , 1 ] be a finite set of rational numbers a/q (here gcd (a,q) = 1) with denominators q bounded 
by (JM / 2 /T 2 • Suppose that 

(7) Bad(W) n [0,M) £ (J -M + ^ M ^ 16 

a/ q z e\Q q q q 
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Then every odd integer K\ < N < K 2 is the sum of seven non-negative cubes. 

Proof. Let N be an odd integer satisfying K\ < N < K 2 - It follows from assumptions (i)-(iii) that 
M = lcrn(W) is squarefree and divisible only by 3 and primes = 5 (mod 6 ). Thus there exists a 
unique integer T £ [0, M) such that 

( 8 ) N = 8T 3 (mod M). 

Suppose N is not the sum of seven cubes. Then, by Lemma 2.1, for each m £ W, we have 
R(T,m) £ [0,?n) \ \e m ■ m,5 m ■ m). Thus T £ Bad(W) n [0 ,M). By (7) there is some rational 
a/q £ 6 such that 




q 


q 


or equivalently 


M , M 

< 8(<jT-aM) 3 < — . 


Moreover, the denominator q is bounded by AL/ 2 K 2 and so 


q 3 N < 


MN M 


2 Ko 


as N < K ‘2 . Hence 


\q 3 N-8(qT-aM) 3 \<M. 


It follows 

□ 


However, by ( 8 ) we have q 3 N - 8(qT - aM) 3 = 0 (mod M). Thus q 3 N = 8(aT - aM ) 
that N is a perfect cube, and so is certainly the sum of seven non-negative cubes. 

We shall mostly apply Proposition 2.2 with the parameter choices given by the following lemma. 

Lemma 2.3. Let K > 10 5 . Let K\ = 9A'/10, K 2 = K. Let 

(9) 


263 -3 292 , /3 

-R l/3 <m< —— A i/d . 


100 100 

Then conditions (iv)-(vi) are satisfied with e m = 0 and S m = 1/10. 


3. Plan for the paper 

The rest of the paper is devoted to understanding and computing the intersections Bad(W’) n 
[0 ,M) appearing in Proposition 2.2. Section 4 collects various properties of remainders and bad 
sets that are used throughout. Section 5 provides justification, under a plausible assumption, that 
the intersection Bad(W) n [0,M) should be decomposable as in (7). Section 6 gives an algorithm 
(Algorithm 1) which takes as input a finite set of positive integers W and an interval [A,B) and 
returns the intersection Bad(>V) n [A,B). We also give a heuristic analysis of the algorithm and 
its running time. Section 7 introduces the concept of a ‘tower’, which a sequence 

(10) W 0 £ Wi £ W 2 £ ••• £ W r = W. 

Letting M,; = lcm(W,:), we prove the recursive formula for computing Bad(Wj) n [0 ,Mf) in terms 
of Bad(Wi-i) n [0, M,_i). This recursive formula together with Algorithm 1 is the basis for a much 
more efficient algorithm (Algorithm 2) for computing Bad(W) n [0, M) given in Section 7. 

In Section 8 we let M* be the product of all primes p < 167 that are = 5 (mod 6 ), and 

(11) W* = { m | M* : 265 x 10 9 < to < 290 x 10 9 } . 

We use a tower and Algorithm 2 to compute Bad(W*) n [0, M*). The actual computation con¬ 
sumed about 18,300 hours of CPU time. 

Section 9 is devoted to proving Theorem 1 for N > (9/10 ) 3 " 8 • J(T « 4.76 x 10 44 , where 
JU = exp(524). The approach is to divide the interval (9/10) 3998 JU < N < into subinter¬ 
vals (9/10)" +1 JU < N < (9/10)"JU with 0 < n < 3997, and apply Proposition 2.2 and Lemma 2.3 
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to prove that all odd integers in the interval (9/10) n+1 J£ r < N < (9/10 ) n J(L are sums of seven 
non-negative cubes. Indeed, we show that given 0 < n < 3997, there is some suitable positive k 
such that the elements of Wo = k • W* satisfy conditions (i)-(iii) (with K\ = (9/10) n+1 J^ and 
K 2 = (9/10)"J^) and that moreover, Bad(Wo) = ft Bad(W*). Thus the results of the huge compu¬ 
tation of Section 8 are recycled 3998 times; on top of this Wo we construct a tower and continue 
until we have found a set W that satisfies the hypotheses of Proposition 2.2, thereby proving 
Theorem 1 for ./V > (9/10) 3 " 8 • J(L. The CPU time for the computations described in Section 9 
was around 10,000 hours. 

The proof of Theorem 1 is completed in Section 10 where a modified strategy is needed to 
handle the ‘small’ ranges (9/10) n+1 Xf < N < (9/10) n J^ with 3998 < n < 4226. Although these 
intervals are small (and few) compared to those handled in Section 9, we are unable to recycle the 
computation of Section 8. This makes the computations far less efficient, though still practical. 
The CPU time for the computations described in Section 10 was around 2,750 hours. 

4. Some Properties of Remainders and Bad Sets 
L emma 4.1. Let m and n be positive integers with k \ m. Then for any real x we have 

q(-, —) = Q (x,m), R (-, —) = - R (x,m). 

V ft ft 1 \ ft ft ) ft 

Let k be a positive integer. For a set Iclwe denote nX - {kx : x e X}. 

Lemma 4.2. Let m and k be positive integers. Let 0 < e < <5 < 1 be real numbers. Then 

Bad(ftm, e,6) = k- Bad(m, £,5). 

Let W be a set of positive integers and for m e W let 0 < £ m < S m < 1 be real numbers. Let 
W — hv * W , £ = (£m)meW ? d — (dm)meW ? § ~ (^m/f:)meW' d — (S m / K )meW' • 

Then 

Bad(W', s', S') = k- Bad( W, e, 5). 

Proof. By (5) and Lemma 4.1, 

x e Bad(«:m, e, 5) R(t, nm) e [0, nm) \ [e • /era, S ■ nm) 

— R(a;, nm) e [0, rn) \ [e • to, 5 ■ m) 

K 

<=> R (x/k, m) e [0, m) \ [e • m, S ■ to) 
x/k e Bad(TO, e, 6) 
x e k ■ Bad(m,£, S). 

This proves the first part of the lemma. The second part now follows from (6). □ 

Lemma 4.3. Given positive integers M\ \ M 2 , we define the ‘natural’ map 
ttm 2 ,Mi : [0,M 2 )—* [0, Mi), x R(x,Mi). 

Then 7Tm 2 ,Mi is surjective, and for any T £ [0, Mi), 

(M 2 /Mi)-1 

n M 2 ,M 1 (T) = U (fc • Mi + T). 
k=0 

Lemma 4.4. Let W2 be sets of positive integers with Wi ^ W2. Write U = W2 - Wi- Let 
LL*i — lcm(Wi). Write 7r — Then 7r(Bad(W2)) £ Bad(Wi) and 

Bad(W 2 ) n [0,M 2 ) = 7r' 1 (Bad(Wi) n [0,Mi)) n Bad (U). 
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Proof. Let if I and let y = n(x) = R(x,Mi). If m e W\ then R (y,m) = R (x,m), as m \ Mi. 
Observe that 

x € Bad(>V2) <=> R(x, to) t \e m ■ to, S m ■ m) for all to e W2 
=► R(x, to ) £ [s m • to, • to) for all to e Wi 
•<=>• R (y, m) i [s m ■ to, S m ■ to) for all to. e Wi 
<=► 2/ € Bad(Wi) ■ 

This shows that 7 r(Bad(W , 2 )) £ Bad(Wi )■ The rest of the lemma easily follows. □ 


5. Gaps and Ripples 


We will soon give an algorithm for computing the intersection 

Bad(W) n [0 ,M) = (n m£W Bad(m)) n [0,M) (M = lcm(W)) 


given a set W that satisfies the conditions of Proposition 2.2. The statement of Proposition 2.2 
(notably (7)) suggests that we are expecting this intersection to be concentrated in small intervals 
around aM/q for certain a/q with relatively small denominators q. In this section we provide an 
explanation for this. The situation is easier to analyze if we make choices of parameters as in 
Lemma 2.3. Thus for this section we fix the choices e m = 0, S. m = 1/10, and hence Bad(m) = 
Bad(?n,0,1/10). We suppose that the elements to e W belong to an interval of the form 


( 12 ) 


263 

100 ' 


L < to < 


292 
loo ’ 


for some L > 0 ( c.f. Lemma 2.3). In fact, we show that if q is large, and if the residues of the 
integers ciM/m are regularly distributed modulo q (in a sense that will be made precise), then 
the intersection Bad(>V) n [0 ,M) contains no points in a certain explicitly given neighbourhood 
of aM/q. Likewise we show for certain a/q with q small, that Bad(>V) n [0, M ) does contain some 
points near aM/q. We stress that the material in this section does not form part of our proof 
of Theorem 1. It does however explain the results of our computations that do form part of the 
proof of Theorem 1, and it lends credibility to them. 


We fix the following notation throughout this section. 

• L is a positive real number; 

• W is a non-empty set of positive integers that belong to the interval ( 12 ); 

• M = lcm(W). 


5.1. Ripples. 


Proposition 5.1. Suppose M > 2000L. Let a/q e [0,1) be a fraction in simplest form with 

1 < q < 9 and 0 < a < q - 1. For 0 < k < 9 - q let 

, 10 , . 292 (k 1\ T 263 (fc+1) 

V ^ 100 V <7 10/ 100 q 

Then ifk <^k and 

(14) M (-M + if k -L, -M + 'I',,, • c Bad(W) n [0, M). 

k =0 V 9 9 > 

This recipe gives 103 disjoint intervals contained in Bad/W) n [0, M) of total length f ■ L where 

(1 d.p.). 


, 261707 

f = _ * 24.9 


10500 


We shall informally refer to the union of intervals (14) as ripple emanating from aM/q 
in the positive direction. The reader will easily modify the proof below to show, under similar 
hypotheses, that there are ripples emanating from the aM/q in the negative direction. 
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Proof. It is easy to check that ipk < ^k for q < 9 and 0 < k < 9 - q. The assumption M > 2000L 
ensures that the 103 intervals are contained in [0,M) and are disjoint, so it is enough to show 
that the intervals are contained in Bad(W). Let a be a real number belonging to the interval 
ipk ' L < a < ■ L. We would like to show that aM/q + a e Bad(?n) for all ?n e W. Let m e W. It 

follows from (12) and (13) that 


(15) 


k 1 \ (k+ 1 ) 

- + — m < 4>k ■ L < a < 'Lfc • L <-to. 

\q 10 / q 


As to | M we can write aM = um with u e Z. Now u = bq + s where 0 < s < q — 1. Thus 


a s 

—M = bm h—to. 

q q 


From (15), 


, (k + s) m a „ _ , (fc+s+ 1 ) 

bm + -—to + — < —M + a < bm + - -m. 

q 10 q q 


Let k + s = qt+v where 0 < v < q - 1. Hence 


(b + t)m + (- + —) m < —M + a < (b + t)m + 
\q 10 / q 


(v+ 1 ) 


TO . 


Observe that 


1 v 1 v +1 

- < — -f - < - 

10 9 10 q 


< 1 , 


as q < 9 and 0 < v < q - 1. Thus Q(aM/q + a, to) -b+t and 


10 


< R 


(aM \ 

^- \-a,mj < to. 


This shows that aM/q+ a e Bad(TO) as required. 


□ 


In the above proposition we showed the existence of ripples emanating from aM/q for q < 9. 
There can also be ripples emanating for aM/q for larger values of q if the sequence of residues 
aM/m in Z/ qL contains large gaps as illustrated by the following proposition. 

Proposition 5.2. Let a/q e (0,1) be a rational in simplest form with q > 11 and 1 < a < q- 1. Let 
(q - 10)/10 < d < q- 1 be an integer, and let s be a non-negative integer satisfying 

263 

(16) s<q-d- 1 , s < ^j( 10 ^ + 10 - q)- 
Suppose 

(17) s+ 1 , s + 2,... ,s + d i |aM/m : m e wj £ Z/gZ. 

Let 

292 s ^ 263 ({s + d+ 1) 1 \ 

100 q 100 \ q 10 / 

Then n < n and 

(-M-H-L,-M-tt-l) £ Bad(W). 

\q q) 

Proof. Let to e >V, and recall that to | M. Thus aM/m is an integer, and hence so is R (aM/m, q). 
By assumption (17), 

R (aM/m, q) * s + 1 , s + 2 ,..., s + d. 

Thus R (aM/m, q) ( (s , s + d + 1). By Lemma 4.1, 

R (aM/q,m) = R(aM,qm) ■ 1/q = R (aM/m,q) ■ m/q. 
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Thus 

(18) 


/ , x / s (s + d+ 1 ) 

K (aM/q,m) f I —m , - m 

\Q Q 


) 


The condition d> (q- 10)/10 implies that 


s (s + d + 1 ) 1 

—to 


Let a belong to the interval 
(19) 


s I (s + d + 1) 

-m < a < 1 

q \ q 


/(£+d+l)_J_\ 

l q 10/ 


We claim that 

R (aM/q- a,m) £ [0,m/10). 

Suppose otherwise: then we can write 

—M - a = bin + r 

q 

where 0 < r < m/10. Thus 

, s , a r , m , (s + d+ 1 ) 

bm + —m < bm + a < —M < bm + a h-< bm + -m 

q q 10 q 

as a satisfies (19). This contradicts (18), and establishes our claim. In fact we have shown that if 
a belongs to the interval (19), then aM/q- a e Bad(?n). 


Suppose now that a belongs to the interval tt • L < a < U ■ L (the second inequality in (16) 
ensures tt < II). To prove the proposition, all we have to show is that a satisfies the inequalities in 
(19) for all m 6 W. However, these follow straightforwardly from the fact that all meW belong 
to the interval ( 12 ). □ 


A few remarks are in order concerning Proposition 5.2 and its proof. 

• For simplicity we have only constructed the first interval in a ripple emanating from aM/q in the 
negative direction. If inequalities (16) are satisfied with a significant margin, then it is possible 
to construct more intervals belonging to this ripple. Likewise, with a suitable modification of 
the assumptions one can also construct a ripple in the positive direction. 

• The first inequality in (16) is imposed merely for simplicity; if it does not hold one can also 
construct ripples emanating from aM/q after suitably modifying the second inequality in (16). 

• The one indispensable assumption in Proposition 5.2 is the existence of a sequence 

s+l,s + 2, ...,s + d 

of consecutive residues belonging to ( Z/qZ ) \ {aM/m : m 6 W) of length d that is roughly 
larger than qj 10. We shall show below that if there is no such sequence, then Bad(>V) contains 
no elements in a neighbourhood of aM/q. 


5.2. Gaps. Let a/q e [0,1] be a rational in simplest form, and let 

&a/q '■ W-^Z/qZ, mn-^a/M/m). 

In view of the above, define the defect d(W,a/q) of W with respect to a/q as the length of 
the longest sequence s + 1, s + 2 ,... ,s + d belonging to (Z/qZ) \ <h a / g (W). As W t 0 , we have 
d{W,a/q) < q. For example, if $ a /q is surjective then d(W,a/q ) = 0, and if $ a / g (W) = (Z/qZ)* 
then d{W,a/q ) = 1 . 
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Lemma 5.3. With notation as above, let cL = d(W,a/q). Let x € R. Then there is some element 
mtW and an integer k such that 



qm 


d + 1 

^r 


Proof. Let u e Z satisfy |w - qx\ < 1/2. We first suppose that d is even. Consider the sequence 

u - d/2, u - d/2 + 1, u- d/2 + 2,..., u + d/2 


of d + 1 elements of Z/qZ. By definition of d, one of these equals 4> a / g (m) for some m e W. Thus 
there is some integer k such that 


aM , 

u - kq 

m 



As \u - qx | < 1 / 2 , the result follows. 

Now suppose that d is odd and qx > u (the case qx < u is similar). Consider the sequence 


u- (d- l)/2, u-(d- l)/2+ 1 , u-(d- l)/2 + 2 ,..., u+ (d+ l)/2 


which again has d + 1 elements, and so there is some rn e W and some integer k such that 


u - 


(d- 1 ) 
2 


aM (d+1) 

< - + kq < u+ -- 

m 2 


Since 0 < qx - u < 1/2, the lemma follows. 


□ 


Lemma 5.4. Let 


Then for all m e W, 


* 


m 


38398 
13875 ’ 


L 

m 


_L_ 725 

m* ~ 38398 


Proof. By (12), the quantity L/m belongs to the interval [100/292,100/263]. We have chosen m* 
so that L/m* is the mid-point of the interval. The lemma follows as 725/38398 is half the length 
of the interval. □ 


Proposition 5.5. With notation as above, let d = d(W,a/q ) and suppose that d < (q- 10)/10. Let 


( 20 ) 

Then 


38398 / 1 (d+l)\ 

/i_ 725 \ 20 2q ) 

( —M - p.L , —M + p,l] n Bad(W) = 0 . 
\q q ) 


A few words are perhaps appropriate to help the reader appreciate the content of the proposi¬ 
tion. We shall suppose that q > 11. If is large compared to q, then we expect that <1 > Qj / 9 is 
close to being surjective and which forces d to be small. If that is the case then qi should be close 
to 38398/(725 x 20) » 2.64. Suppose now that #W is large, but that q is much larger. Suppose 
also that the residues in the image 4 > a / g (W) are ‘randomly’ distributed in Z/qZ. The quantity d 
measures how large the gaps between these residues in the image can be, and we expect that d 
should be around q/ffW. We therefore expect that /./, « (38398/725) (^ - 2 ._^ vv ). We see that /./, 
should be positive if W has much more than 10 elements. 


Proof of Proposition 5.5. The assumption d < (q - 10)/10 ensures that /i is positive. Let y e 
( aM/q-nL,aM/q+p,L ). We would to like to show that there is some m e W such that y i Bad(m). 
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Write y = aM/q + (3 where |/3| < /. iL. Letting x = 1/20 - (3/m* in Lemma 5.3, we deduce the 
existence of some integer k and some element raeW such that 


(3 aM , 1 

— +- + k - 

m* qm 20 


(d+ 1 ) 
2 q 


Thus 


(3 aM 

— +-+ k - 

m qm 


1 

20 


(d + 1) 

< i-- + 

- 2 q 


]3_ 

m* 


1 

TO 


Using \(3\ < pL, Lemma 5.4 and the definition of p in (20), we see that 


(3 aM , 

— +-+ k - 

m qm 


1 

20 


20 


Thus y = aM/q + (3 belongs to the interval -km+ (0,to/ 10), showing that y £ Bad(m) as required. 

□ 


6. A First Approach to Computing Bad(W) 

In this section W is a finite set of positive integers m. Associated to each to e W are real 
numbers 0 < £ m < S m < 1. We shall write e - (Sm/meU and § = (<5m)meW- 


Lemma 6.1. Let A < B be real numbers. For m e W, Zet 

<? TO = Q(A,m), r m =R(A,rn). 

/a/ Suppose r n e[e n -n,8 n -ri) for some nt W. Write A' = min((< 7 „ + 5„) • n, B). Then 

Bad(W) n [A, B) = Bad(W) n [A', B). 

(b) Suppose r m [e m ■ to, S m ■ to) for all m e W. Define 


( 21 ) 


A 


m 


Then 


(q m + £ m )-m ifr m <e m ■ to 
(< 7m + 1 + e m ) • TO ifr m >Sm ■ m, 


A' = min (B,min( A m ) m€V v). 


Bad(W)n[A,B) = (Bad(W) n [A',B)) u [A, A'). 


Proof. Suppose n e W satisfies r n € [e n • n, S n -n), and let A' be as in (a). By (5) we have 

( q n ■ n + [s n ■ n, 6 n ■ n)) n Bad(n) = 0 . 

Observe that [A, A') Qq n -n+ [s n -n,5 n - n ) and [A, A') £ [A, _B). Part (a) follows. 

Suppose now that r m £ [e m • rn,S m ■ to) for all to e W, and let A' be as in (b). It is easy to 
check that R(A",to) i [e m • rn. S rn - to) for all A" e [A,A m ). From this we see that [A, A') £ 
n m€ wBad(TO,e m ,(I m ) =Bad(W’). Part (b) follows. □ 


Lemma 6.1 immediately leads us to the following algorithm. 

Algorithm 1. To compute Bad(W) n [ A,B ) as a disjoint union of intervals U itj? I- 
Input: A, B, W, e, 8. 

Initialize J? *- 0 . 

Repeat the following steps until A = B. 

(a) Loop through the elements to 6 W computing q m = Q(A,to), r m = R(A,to). 

(b) If there is some n e W such that e n - n <r n < S n ■ n then 

A «- min ((q n + S n ) ■ n, B ) 

and go back to (a). 

(c) Otherwise, let A' be as in (21). Let JA *- J^u{[A,A')} and then A <- A'. Go back to (a). 

Output: J?. 
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A heuristic analysis of Algorithm 1 and its running time. Let x e [0, M) and recall that 
R(x,r?z) e [0 ,to). Moreover, x e Bad(m) if and only if R (x,m) e [0 ,to) \ [e rn ■ m,S m ■ m). Thus 
the ‘probability’ that x belongs to Bad(m) is 1 - (S m ~ £ m )- Assuming ‘independence of events’ 
we expect that the total length of intervals produced by Algorithm 1 is 

(22) (B-A)- ]1 + 

meW 

To analyse the running time, we shall suppose parameter choices as in Lemma 2.3: namely 
e m = 0 and 6 m = 1/10 for all m e W. Moreover, we shall suppose that the elements of m e W 
belong to an interval (12) for some large positive L. By the above, the expected total length 
of the intervals produced by Algorithm 1 is (B - A) ■ 0.9# w . Moreover, we suppose that W is 
sufficiently large so that the length of the output should be negligible compared to B-A; this 
should mean that step (c) is relatively rare. We will estimate the expected number of times we 
loop through steps (a), (b). Note that in step (b), A is increased by 0.1 -n-r n . The remainder 
r n = R(A, n) belongs to [0,0.1 • n). We regard the increase as a product (0.1 - r n /n) ■ n. Treating 
r„/n as a random variable uniformly distributed in [ 0 , 0 . 1 ) and n as a random variable uniformly 
distributed in interval (12), we see that the expected increase is 0.05-(2.63 + 2.92)A/2 = 0.13875-A. 
A standard probability theory argument that we omit tells us that the expected number of times 
the algorithm loops through steps (a), (b) is roughly 

(B - A)/(0.13875A) » 7 (B-A)/L. 

We now suppose that AT is very large, and we would like to compute the intersection Bad(W) n 
[0, M) for some set W where we hope that the hypotheses of Proposition 2.2 and Lemma 2.3 are 
satisfied. In particular, we take L = A' 1 / 3 . The number of steps should be around 7M/A' 1 / 3 . We 
have to choose W so that M = lcm(>V) is much larger than K (see (7) and just above it). Thus the 
number of steps to compute Bad(W’) is much greater than A' 2 / 3 . For K = exp(524), the expected 
number of steps is larger than 10 15 °, which makes the computation entirely impractical. 

7. A Refined Approach to Computing Bad(W): The Tower 

In this section we let W be a set of positive integers with M = lcm(W). Let Mo, Mi, M 2 , ■■■, M r 
be positive integers such that Mj | M, + i and M r = M. Write pi = Mj+i/Mj. In our later computa¬ 
tions the pi will be primes, but we need not assume that yet. Let 

Wi = {m e W : m \ Mi}. 

We suppose that M, : = lcm(W,;). Write Ui = W /+1 \ W,. Recall (Lemmas 4.3 and 4.4) that we have 
natural surjections 7 ^ : [0 ,Mj) -*• [0,M, ; ) whenever j > i, and that these restrict to give maps 
(not necessarily surjections) Bad(Wj) Bad(>Vj). We shall refer to the sequence of inclusions 
(10) as a tower leading up to Bad(>V), and use this to compute Bad(>V). 

Lemma 7.1. Let 0 < i < r - 1. Suppose J’i is a finite set of disjoint subintervals of [0,M,) such 
that 

Bad(Wi) n [ 0 , Mf) = \J I. 


Bad(W i+ i) n [0, M i+1 ) = (J \J ((k ■ + I) n Bad {U}) . 

k=0 

Proof. This is immediate from Lemmas 4.3 and 4.4. □ 

Lemma 7.1 immediately leads us to the following algorithm. 
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Algorithm 2. The following computes a finite set JA - jA r of subintervals of [0, M) such that 

Bad(W)n[0,M) = U us I. 

Input: Wo, • ■ •, W r = W, e, 5. 

Initialize: J 5 o to be the set of disjoint intervals whose union equals Bad(Wo) n [0,Mo), and which 
is computed using Algorithm 1. 

Initialize: i *- 0. 

Repeat the following steps until i = r. 

(a) A + i 0 . 

(b) for I 6 jAi and k e {0,... ,Pi~ 1}, compute, using Algorithm 1, a finite set JA' of subintervals 

of [0, Mi+ 1 ) such that ( k ■ + I) n Bad (Uf) - U i'^' I'; let J^+i *- J^+i u y. 

(c) i <- i+ 1. 

Output: JA = J* r . 

A heuristic analysis of Algorithm 2 and its running time. We shall suppose, as in 
Lemma 2.3, that e m - 0 and 6 m = 1/10 for all m € W. Write rq = #W,;. We assume that 
the elements of W t , W; belong to an interval of the form [263L/100,292L/100] for some large L. 
By our previous analysis, we expect that we can compute JA$ in roughly 7M 0 /L steps. The total 
length t(jA 0 ) of the intervals in J’e should roughly be 0.9"°Mo. In Step (b) of the algorithm, we 
will replace each / 6 J^o with po intervals of the same length, and then apply Algorithm 1 to each. 
Thus we expect that the number of steps to compute ,JA\ to be roughly 

7p 0 • 0.9"° • M 0 7Mi • 0.9 n ° 

L ~ L ' 

The total length of the intervals in should be roughly Mi ■ 0.9" 1 . It is now apparent that the 
total number of steps should be around 

(7/L) ■ (M 0 + Mi • 0.9”° + M 2 • 0.9 ni + - + M r ■ 0.9 n ’- 1 ). 

8 . A Large Computation 

Let M* be the product of all primes p < 167 that are = 5 (mod 6 ), and W* is as in (11). In 
this section we compute Bad(W*) n [0,M*), using a tower and Algorithm 2. As explained in the 
plan (Section 3), the result of this computation will be reused again and again in Section 9. Let 

M 0 = 5 x 11 x 17 x 23 x 29 x 41 x 47 x 53 x 59 x 71 x 83 x 89, 

which is the product of the primes < 100 that are s 5 (mod 6 ). Let 

Mi = 101 • M 0 , M 2 = 107 • Mi, M 3 = 113 • M 2 , M 4 = 131 • M 3 , 

M 5 = 137 • M 4 , M 6 = 149 • M 5 , M* = M 7 = 167 • M 6 . 

We let 

W, = { m | Mi : 265 x 10 9 <m< 290 x 10 9 }. 

Thus Wo £ ••• £ W 7 = W*. We checked that M,; = lcm(Wj). Table 1 gives the cardinalities of the 
Wj. We use this tower and Algorithm 2 to compute Bad(W’ 1 ') n [0,M*). By our heuristic in the 
previous section, the number of steps needed for this computation should very roughly be equal 
to 6.0 x 10, which is the sum of the entries of the table’s third column. It appears from this 
estimate that the computation can be done in reasonable time. 

We wrote simple implementations of Algorithms 1 and 2 for the computer algebra system Magma 
(Bosnia et ah, 1997). We divided the interval [0,M 0 ) into 59,000 subintervals of equal length 
and ran our program on each of these intervals \Ak-i,Ak) successively computing Bad(W ? ) n 
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Table 1. The M* and the Wi are given at the beginning of Section 8. The third 
column gives an estimate for the number of steps needed to compute Bad(Wi) n 
[0, Mi) from Bad(Wi-i) n [0, M,_i) according to the heuristic analysis at the end 
of Section 7. 


i 

logio (Mi) (! d.p.) 

ni = #W 4 

7Mi ^ 1 ( 2 s-f.) 

0 

18.3 

16 

1.4 x 10 9 

1 

20.3 

38 

2.6 x 10 9 

2 

22.3 

83 

2.7 x 10 10 

3 

24.4 

149 

2.7 x 10 10 

4 

26.5 

250 

3.4 x 10 9 

5 

28.6 

401 

1.1 x 10 7 

6 

30.8 

620 

2.0 x 10 2 

7 

33.0 

911 

3.2 x 10" 6 


7 Tj- \([Ak-i,A}.)) for i = 0,...,7. Our computation was distributed over 59 processors (on a 64 core 
machine with 2500MHz AMD Opteron Processors). Note that 

59,000 

Bad(W’ i ) n [0, Mf) = \J Bad(W;) mr^A^, A k ))- 

k= 1 

thus our computation gives us a decomposition of Bad(>Vi)n[0, Mj) as a union of disjoint intervals. 
The total CPU time for the computation was around 18,300 hours, but as we distributed the 
computation over 59 processors, it was over in less than two weeks. 

Lemma 8.1. There are sequences (Bj)^f, (Cj)j=i contained in [0 ,M*\ such that 

B\ < Ci < £?2 < Cl < ••• < i?854 < Cs54 

and 

854 

Bad(VW)n[0 ,M*) = (J [B^C,), 

3= 1 

854 a 

with total length Y,( C j ~ D .i) = 20382195221000—. 

j =i 10 


Proof. As indicated by Table 2, our computation gives Bad(W*) n [0 ,M*) as a union of 861 
intervals disjoint subintervals of [0 ,M*). Among these there are 7 pairs of the form [a,/3)u[/3,7), 
where the values of fd are of the form /3' • M */59000 with 

P' = 7375, 14750, 22125, 29500, 36875, 44250, 51625. 

These subdivisions are clearly a result of our original subdivision of interval [0, Mq ) into 59000 
subintervals of equal length. We simply replace the pairs [«,/?) u [/3,y) with [a, 7) so that 
Bad(yV’’ 1 ') n [0 ,M*) is expressed as a union of 854 intervals. This simplification of course pre¬ 
serves the total length of intervals. □ 


8.1. Remarks and Sanity Checks. Our computations are done with exact arithmetic. The 
reader will note by looking back at Algorithms 1 and 2 (and recalling that all e m = 0 and S m = m/10) 
that the end points of the intervals encountered will be rationals with denominators that are 
divisors of 10, except for the Ak appearing in our original subdivision which have denominators 
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Table 2. Some details for the computation described Section 8. The second 
column gives where is a disjoint collection of intervals such uj^ = 

Bad(yV’ J ) n [0, Mj). The third column gives the total length £, of these inter¬ 
vals. The fourth column gives (2 s.f.) the ratio A/M,;. According to the heuristic 
at the end of Section 6, this ratio should approximately equal 0.9"* which is 
given (2 s.f.) in the last column (here n,; = #W,: as in Table 1). We explain the 
discrepancy between the last two columns in Subsection 8.1. 


i 


li = f(Bad(Wj) n [0, Mi)) 

iil Mi 

0.9"* 

0 

23,458,002 

365,300,497,739,376,385^ 

1.85 x 10' 1 

1.85 x 10' 1 

1 

553,209,618 

3,625,384,986,862,035,664^ 

1.82 x 10 -2 

1.82 x 10 -2 

2 

1,106,375,245 

3,313,998,145,602,553,709^ 

1.56 x 10" 4 

1.59 x 10” 4 

3 

209,982,392 

350,826,426,611,537,217^ 

1.46 x 10" 7 

1.52 x 10" 7 

4 

1,062,201 

1,076,402,154,947,217^ 

3.41 x 10“ 12 

3.64 x 10“ 12 

5 

904 

20,663,973,893,432^ 

4.78 x 10~ 16 

4.48 x 10“ 19 

6 

870 

20,504,346,087,851^ 

3.19 x KT 18 

4.27 x 10“ 29 

7 

861 

20,382,195,221, OOOyjy 

1.90 x 10" 20 

2.07 x 10“ 42 


that are divisors of 59,000. As a check on our computations, we verify that our results for 
Bad(W*) n [0, M*) are consistent with Proposition 5.1. The set W* satisfies 

min(VW) = 265,024,970,473 max(W*) = 289,916,573,827. 


We take L = min(>V*) • 100/263. It turns out that L > max(W*) • 100/292. Thus W* is contained 
in the interval (12) for this value of L. Proposition 5.1 yields a total 103 intervals of the form 
( aM*/q + ipk • L, aM*/q + T*, • L ) that must be contained in Bad(>V*) n [0, M*). We checked that 
each of these is contained in one of the 854 intervals produced by our computation. It is instructive 
to compare the fourth and fifth columns of Table 2. According to our heuristic, the total length 
f'(Bad(yV’;)n[0, Mi)) should be around Mj-0.9"* (with n,; = #>Vj) and therefore we expect the two 
columns to be roughly the same. From the table, we see that this heuristic is remarkably accurate 
for 0 < i < 4, and extremely inaccurate for i > 5. An explanation for this is provided by the ripples. 
The total length of the intervals contained in Bad(W,) n [0, Mj) produced by Proposition 5.1 is 
« 24.9 L. Now 


24.9L 

M s 


= 5.8 x 10“ 17 , 


24.9L 

M e 


= 4.0 x 10 19 , 


24.9L 

M r 


= 2.3 x 10" 21 


(1 s.f.), 


which does provide an explanation for the discrepancy between the too columns. Proposition 5.2 
(with W* and M* in place of W and M) produces 172 intervals with 11 < q < 100 with total length 
« 17.8 L. We checked that each of these is also contained in one of the 854 intervals produced by 
our computation. 


According to the overall philosophy of Section 5, the set BadClV*) n [0,M*) should be con¬ 
centrated in short intervals around rational multiples ( a/q ) • M* with q small. To test this, we 
computed, using continued fractions, the best rational approximation to (Bi + Q)/(2M*) with 
denominator at most 10 20 , for 1 < i < 854. The largest denominator we found was 42. 
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The reader is probably wondering, given that we are employing 59 processors, why we have 
subdivided [0,Mo) into 59,000 intervals instead of 59 intervals. This was done purely for memory 
management reasons. A glance at Table 2 will show the reader that there is an explosion of 
intervals at levels i = 1, 2, 3. By dividing [0, Mo) into 59,000 subintervals, we only need to store 
roughly 1/59000-th of the intervals appearing at levels i at any one time per processor, and so 
only need to store around 1/1000-th of these intervals in the memory at any one time. 

9. Proof of Theorem 1 for N > (9/10) 3 " 8 • exp(524) « 4.76 x 10 44 

The reader might at this point find it helpful to review the first paragraph of Section 2 as well as 
the plan in Section 3. Let X = exp(524). In this section we prove Theorem 1 for N > (9/10 ) 3998 X. 
We shall divide the interval (9/10 ) 3998 X <N <X into subintervals (9/10)" +1 JT < N < (9/10)"JT 
with 0 < n < 3997. We apply Proposition 2.2 and Lemma 2.3 to prove that all odd integers in the 
interval (9/10)" +1 J£T < N < (9/10)"J£T are sums of seven non-negative cubes. 

Lemma 9.1. Let 0 < n < 3997. Let K = (9/10)" • X. There exists an integer k that satisfies 

(a) k is squarefree; 

(b) 3 | k 

(c) k/ 3 is divisible only by primes q = 5 (mod 6) that satisfy q > 167; 

(d) k belongs to the interval 

263 K 1 ! 3 ^ 292 K 1 ' 3 

( ' 23 ' 265 ’ 10 11 “ h ~ 290 ’ 10 11 ' 

Proof. We proved the lemma using a Magma script. Let I±, I 2 be the lower and upper bounds for 
k in (23). If I 2 < 10 7 then our script uses brute enumeration of integers in the interval \I\,I 2 \ to 
find a suitable k. Otherwise, the script takes r to be a product of consecutive primes = 5 (mod 6) 
starting with 173 up to a certain bound, and keeps increasing the bound until / 2 /r < 10'. It then 
loops through the integers /i/3r < p, < J 2 /3t until it finds one such that n = 3 pr satisfies conditions 
(a), (b), (c). □ 

Remark. For n = 3998, the interval in (23) is 7481.6 •••<«< 7590.5 ..., which is too short for the 
existence a suitable k. This is also the case for most values of n that are > 3998. 

Lemma 9.2. Let 0 < n < 3997 and let n be as in Lemma 9.1. Let W* and M* be as in Lemma 8.1. 
Let 

Wo = {k- m* : m* e W*}, M 0 = lcm(W 0 ) = kM*. 

Let £m = 0 and S m = 1/10 for all m e Wo- Then in e Wo satisfy the conditions (i)-(vi) of Section 2, 
where 

Ki = (9/10)" +1 • JT, K 2 = (9/10)" • X. 

Moreover, 

854 

(24) Bad(W’o) n [0, M 0 ) = \J [k ■ Bj , k • Cj), 

3 =1 

where the Bj and Cj are as in Lemma 8.1. 

Proof. All rri* € W* are squarefree and divisible only by primes q < 167 satisfying q = 5 (mod 6). 
Thus conditions (i)-(iii) of Section 2 are satisfied by m e Wo- As we are taking e m = 0 and 
- 1/10, to verify conditions (iv)-(vi) we may apply Lemma 2.3. For this we need only check 
(9) holds for to e Wo, where K = K 2 . This immediately follows from (23) and the fact that 
W* c [265 x 10 9 ,290 x 10 9 ]. 

Finally, by Lemma 4.2, 

Bad(W’o) n [0, M 0 ) = k ■ (Bad(W*) n [0, M*)). 

Lemma 8.1 completes the proof. □ 
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Our Magma script for proving Theorem 1 in the range K\ < N < proceeds as follows. We 
inductively construct a tower Wo c Wi c W 2 c •••. Observe that 

l(Bad(Wo) n [0, Mg)) _ l(Bad(>W) n [0, M*)) ^ x gQ x 1() -20 
M 0 M* 

thus the computation of the previous section has already substantially depleted the interval 
[0,Mo). Given Wi, and Mi, we let pi be the smallest prime = 5 (mod 6) that does not divide M t 
and let Mj+i = jyM>. The script then writes down positive integers m belonging to the interval 
(9), such that to | M i+ 1 and 3 pi \ m. It is not necessary or practical to find all such integers, 
but we content ourselves with finding around 3 log(p,)/log(0.9 _1 ) of them; we explain this choice 
shortly. These to will form the set Ui and we take W , + 1 = W t u Hi. The script then applies our 
implementation of Algorithm 2 to compute Bad(W J+ i) n [0, M i+l ) as a union of disjoint intervals. 
Our heuristic analysis of Algorithm 2 suggests that f(Bad(Wj+i)n[0,Mj + i)) should roughly equal 
p,:-0.9 #Wi T(Bad(Wj)n[0, Mi)). We desire the total length of the intervals to decrease in each step 
of the tower, so we should require > log(pj)/ log(0.9 _1 ). Experimentation suggests that requir¬ 
ing « 3 log (pi)/ log(0.9 _1 ) provides good control of both the total length of Bad(Wi) n [0, Mi) 
and the number of intervals that make it up. Our script continues to build the tower and compute 
successive Bad(Wi) n [0, A/,;) until it finds W = W,: and M - Mi that satisfy (7) for some set of 
rationals 6 c [0,1] with denominators bounded by \J M/2K. Specifically, once Mi > 2K , for each 
of the disjoint intervals [a,f3) that make up Bad(W) n [0, M), the script uses continued fractions 
to compute the best rational approximation a/q to (a + /3)/2M with q < \J M/2K, and then checks 
whether [a,/3) £ ( aM/q- \J M / 16/ q , aM/q+ \J M/!&/q). The script continues constructing the 
tower until this criterion is satisfied for all the intervals making up Bad(W). It then follows 
from Proposition 2.2 that all odd integers in the range ■ (9/10) n+1 < N < ■ (9/10)” are 

sums of seven non-negative cubes. We again distributed the computation among 59 processors on 
the afore-mentioned machine, with each processor handling an appropriate portion of the range 
0 < n < 3997. The script succeeded in finding an appropriate W for all n in this range. The entire 
CPU time was around 10,000 hours, but as the computation was distributed among 59 processors 
the actual time was around 7 days. 


We give more details for the case n = 0. Thus K = = exp(524), and we would like to show, 

using proposition 2.2 that all odd integers 9A'/10 < N < K are sums of seven non-negative cubes. 
The routine described in the proof of Lemma 9.1 gives the following suitable value for k: 


k = 3 x 173 x 179 x 191 x 197 x 227 x 233 x 239 x 251 x 257 x 263 x 269 x 281 x 293 

x 311 x 317 x 347 x 353 x 359 x 383 x 389 x 401 x 419 x 431 x 443 x 207443. 


Table 3 gives some of the details for the computation. We take W = W 48 - Then #>V = #Wo + 
Y, = 9943, and 

^(Bad(W) n [0,M)) =12459371373955496388240157141404031514014113708989680551759 

37887691670913319978- « 1.25 x 10 78 . 

2 

In comparison, 

M = M 48 « 1.64 x 10 235 , K = 3.72 x 10 227 , 

The set © as in (7) turns out be precisely the set of 171 rationals a/q e [0,1] with denominators 
q belonging to 


1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 18, 19, 21, 24, 26, 28, 30, 36, 42. 

As a check on our results, we apply Proposition 5.2 to show that there is an interval close to 
(a/ 42) • M for 1 < a < 41 with gcd(a,42) = 1. Our W and M satisfy the hypotheses of Section 5 
with L = A' 1 / 3 . Note that 3 | to | M for all m e W. As M is squarefree, we have 3 \ (M/m). 
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Table 3. This table gives details for the computation for the case n = 0. For 
each i > 1, our script computes Bad(W;) as a disjoint union of subintervals of 
[0, Mi). The number of intervals is given in the fourth column. The fifth column 
gives, to 3 significant figures, the ratio ti/Mi where li = ^(Bad(Wi) n [0,Mj)). 


i 

Pi-i 

mi-i 

number 

of 

intervals 

lit Mi 

0 

- 

- 

854 

1.90 x 10~ 2 ° 

1 

449 

174 

775 

3.73 x 10- 23 

2 

461 

175 

745 

7.94 x 10“ 26 

3 

467 

176 

740 

1.70 x 10“ 28 

4 

479 

176 

735 

3.54 x 10~ 31 

5 

491 

177 

732 

7.20 x 10~ 34 

6 

503 

178 

730 

1.42 x 10“ 36 

7 

509 

178 

730 

2.80 x 10~ 39 

8 

521 

179 

730 

5.38 x 10^ 42 

9 

557 

181 

730 

9.65 x 10~ 45 

10 

563 

181 

731 

1.71 x 10" 47 

11 

569 

181 

730 

3.01 x 10- 50 

12 

587 

182 

729 

5.13 x 10“ 53 

13 

593 

182 

729 

8.64 x 10" 56 

14 

599 

183 

729 

1.44 x 10~ 58 

15 

617 

183 

729 

2.34 x 10~ 61 

16 

641 

185 

729 

3.64 x 10~ 64 

17 

647 

185 

729 

5.63 x 10^ 67 

18 

653 

185 

729 

8.62 x 10" 7 ° 

19 

659 

185 

729 

1.31 x 10" 72 

20 

677 

186 

729 

1.93 x 10~ 75 

21 

683 

186 

729 

2.83 x 10" 78 

22 

701 

187 

729 

4.04 x 10~ 81 

23 

719 

188 

729 

5.61 x 10^ 84 

24 

743 

189 

729 

7.55 x 10“ 87 


i 

qi-i 

mi-i 

number 

of 

intervals 

til Mi 

25 

761 

189 

729 

9.93 x 10“ 9 ° 

26 

773 

190 

729 

1.28 x 10- 92 

27 

797 

191 

729 

1.61 x 10" 95 

28 

809 

191 

729 

1.99 x 10- 98 

29 

821 

192 

729 

2.43 x 10~ 101 

30 

827 

192 

729 

2.93 x 10” 104 

31 

839 

192 

729 

3.50 x 10- 107 

32 

857 

193 

729 

4.08 x 10" 110 

33 

863 

193 

729 

4.73 x 10” 113 

34 

881 

194 

729 

5.36 x 10” 116 

35 

887 

194 

729 

6.05 x 10~ 119 

36 

911 

195 

729 

6.64 x 10~ 122 

37 

929 

195 

729 

7.14 x 10~ 125 

38 

941 

195 

729 

7.59 x 10” 128 

39 

947 

196 

729 

8.02 x 10~ 131 

40 

953 

196 

729 

8.41 x 10” 134 

41 

971 

196 

729 

8.66 x 10” 137 

42 

977 

197 

729 

8.87 x 10- 140 

43 

983 

197 

729 

9.02 x 10^ 143 

44 

1013 

198 

729 

8.91 x 10” 146 

45 

1019 

198 

729 

8.74 x 10~ 149 

46 

1031 

198 

729 

8.48 x 10” 152 

47 

1049 

199 

729 

8.08 x 10" 155 

48 

1061 

199 

729 

7.62 x 10~ 158 


Moreover, all the prime divisors of M/3 are = 5 (mod 6). It follows that gcd(aAf/m,42) = 1 for 
all m e W. Let q - 42 and s = d = 5 in Proposition 5.2; hypothesis (16) is trivially satisfied. 
Now s + 1,..., s + d are the integers 6, 7, 8, 9, 10 and none of these are coprime to 42. Thus 
condition (17) is also satisfied. By Proposition 5.2, for each 1 < a < 41 with gcd(a,42) = 1 we have 

(25) (-M- ■ K 1/3 , — M - — • A' 1/3 ') £ Bad(W). 

V ’ V 42 10500 42 210 ) V ’ 
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One of the 729 intervals that make up Bad( W) is [it, v ) where the end points u, v are 
u =3895173640423584874713349032421520960246664873653293975537400307 
5224581570150578968661382487115397667257923729694373737120676906 
3932017310777324617938079775100516093231041460322490961793995991 

410145937421686204642056677472293123392066—, 

10 ’ 

v =3895173640423584874713349032421520960246664873653293975537400307 
5224581570150578968661382487115397667257923729694373737120676906 
3932017310777324617938079775101078697615391607190077739469387928 
238665618669912989320140106379011502569660, 


and we checked that the interval in (25) with a = 1 is contained in [it,i/). It is also interesting to 
note how close the two intervals are in length: the ratio of the lengths of the two intervals is 

(4471/10500- 73/210) • K 1 ^ ^ Q ^ 
v - u 

which illustrates how remarkably accurate our Proposition 5.2 is. 

10. Completing the Proof of Theorem 1 

It remains to apply Proposition 2.2 to the intervals (9/10)™ +1 JC < N < (9/10)"JC with 3998 < 
n < 4226. We write K = K 2 = (9/10)™JC and K\ = (9/10) n+1 J#b It is no longer practical to use 
the choices in Lemma 2.3 as the interval in (9) is too short to contain many squarefree in whose 
prime divisors are 3 and small primes = 5 (mod 6). The interval in (9) is a result of imposing 
the uniform choices s m = 0 and S rn = 1/10. Instead we consider integers m satisfying conditions 
(i)- (iii) of Section 2 but belonging to the (much larger) interval 

(26) —KW^mK—K 1 / 3 . 

5 5 

For each such m we take e m = er'/1000 and 8 m = d'/lCHX) where s', S' are integers with s' and S' 
respectively as small and as large as possible such that the conditions (v), (vi) of Section 2 are 
satisfied. We only keep those values of m for which 

(27) 0 < Em < Sm < 1, S m ~ £m — 1/^0; 

an elementary though lengthy analysis in fact shows that the inequalities in (27) together with 
(v) and (vi) force m to belong to the interval (26). Note that the set Bad(m, s m , S m ) has ‘relative 
density’ 1 - S m + £ m in the restriction S m - e m ^ 1/20 ensures that this relative density is not 
too close to 1, and that therefore rn makes a significant contribution to depleting the intervals in 
Algorithms 1 and 2. 

We choose a prime q = 5 (mod 6), depending on K, and let 

M 0 = 3-5-11 . q 

which is the product of 3 and the primes < q that are = 5 (mod 6). Let Wo be the set of positive 
integers dividing Mq and satisfying the above conditions. We found experimentally that for each 
n in the above range it is always possible to choose q so that 

Mo = lcm(W 0 ), n (l-S m +£m)<l/5, log 10 (M 0 /# 1/3 )<7.5. 

meW o 

The inequality FlmeWo(l -< ^m+£m) - 1 /5 indicates that f(Bad(Wo)n[0,M 0 )) should heuristically 
be at most Mo/5 which means that this is a good first step at depleting the interval [0,Mo). 
The other inequality indicates that we can compute Bad( Wo) n [0, Mo) in a reasonable number of 
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steps according to the heuristic following Algorithm 1. We let po be the first prime = 5 (mod 6) 
that is > q, and p\ the next such prime and so on. We let M i+ \ - piA'U and construct a tower 
as before. We stop once Bad(W ? ) n [0. M,) satisfies the criterion of Proposition 2.2. Our Magma 
script succeeded in doing this for all n in the range 3998 < n < 4226. The total CPU time was 
around 2750 hours, but the computation was spread over 59 processors so the actual time was less 
than 2 days. 

We give a few of details for the computation for the value n = 4226. The final M is the product 
of 3 and the primes p = 5 (mod 6) that are < 227. The final W has 8083 elements. It turns out 
that Bad(W) n [0, M) consists of 305 intervals and that ^(Bad(W) n [0« 2.24 x 10~ 32 . 
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